DTP-MSU/06-20 



Cylindrically symmetric solitons in Einstein— Yang— Mills theory 
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Recently new Einstein- Yang-Mills (EYM) soliton solutions were presented which describe super- 
conducting strings with Kasner asymptotics (hep-th/0610183). Here we study the static cylindrically 
symmetric SU(2) EYM system in more detail. The ansatz for the gauge field corresponds to su- 
perposition of the azimuthal B v and the longitudinal B z components of the color magnetic field. 
We derive sum rules relating data on the symmetry axis to asymptotic data and show that generic 
asymptotic structure of regular solutions is Kasner. Solutions starting with vacuum data on the 
axis generically are divergent. Regular solutions correspond to some bifurcation manifold in the 
space of parameters which has the low-energy limiting point corresponding to string solutions in 
flat space (with the divergent total energy) and the high-curvature point where gravity is crucial. 
Some analytical results are presented for the low energy limit, and numerical bifurcation curves are 
constructed in the gravitating case. Depending on the parameters, the solution looks like a straight 
string or a pair of straight and circular strings. The existence of such non-linear superposition of 
two strings becomes possible due to self-interaction terms in the Yang-Mills action which suppress 
contribution of the circular string near the polar axis. 

PACS numbers: 04.20.Jb 



I. INTRODUCTION 

The discovery by Bartnik and McKinnon [| of a 
particle-like solution in the Einstein- Yang-Mills (EYM) 
theory and construction of the corresponding black 
holes [2| stimulated wide research of Einstein- Yang-Mills 
(EYM) solitons and their generalization in supergravities 
and other related models (for a review see [3|]). Between 
them, cylindrically symmetric four-dimensional configu- 
rations remained relatively less explored. Useful general 
analysis can be found in the study of Einstein Maxwell 
geons [1, [H, the Melvin magnetic universe Q and the 
abelian cosmic strings [7(. Other configurations such as 
the Yang-Mills vortex in 4+1 dimensions [§], cylindri- 
cal black holes [§, [l(| and some numerical solution in 
are also relevant. Recently a new class of the EYM 
cylindrically-symmetric solutions was reported [12[ de- 
scribing the Einstein- Yang-Mills strings. The purpose of 
this paper is to study the cylindrically symmetric four- 
dimcnsional EYM system in more detail. 

Our ansatz for the gauge field contains two functions 
which describe the longitudinal and the azimuthal mag- 
netic field components. The stress-energy tensor does 
not exhibit the boost symmetry as in the case of cosmic 
strings, so the metric is parametrized by three different 
functions of the radial variable. We analyze first the so- 
lutions in flat space-time imposing the regularity condi- 
tions on the polar axis. Generically, the integral curves 
starting with such boundary data on the axis diverge at 
finite distance exhibiting the square-root singularity of 
the Yang-Mills field, reminiscent of that in the spheri- 
cally symmetric case [ljj]. There exist, however, the bi- 
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furcation curve in the parameter space which correspond 
to regular solutions. The energy density does not fall 
sufficiently fast at infinity, so the total energy diverges. 

When gravity is switched on, field configurations be- 
come effectively compact and the metric is asymptoti- 
cally Ricci-flat. It is, however, not Minkowskian, but 
Kasner. We derive the integral sum rules relating the 
boundary values of the metric functions on the polar axis 
with the Kasner parameters which show the absence of 
asymptotically flat configurations of non-zero mass. We 
also present a (non-rigorous) proof of existence of regular 
solutions with Kasner asymptotic. Numerical integration 
reveals the following typical behavior of the solutions. In 
the vicinity of the polar axis the space-time is flat. With 
growing distance the longitudinal magnetic field is dom- 
inant within a certain region. The YM field there is ef- 
fectively U(l), while the metric is approximately Melvin. 
With the further growing distance the azimuthal com- 
ponent of the YM magnetic field comes into play and 
the YM non-linearity becomes manifest. In this region 
the solution may be interpreted as describing the circu- 
lar magnetic string. Finally, in the asymptotic region the 
metric acquires the boost symmetry in the £ — </? plane ex- 
actly in the same way as the metric of the cosmic string is 
boost-symmetric in the t — z plane. Moreover, the met- 
ric component g zz asymptotically tends to zero, so the 
space-time becomes effectively 1 + 2 dimensional. This 
effect is similar to the asymptotical contraction of the 
azimuthal dimension in the Melvin solution. Thus, our 
solution may be loosely interpreted as describing the sys- 
tem of interacting straight and circular magnetic strings 
in the four-dimensional EYM theory. 

The paper is organized as follows. In Sec. II we intro- 
duce the gauge field ansatz and discuss the symmetries 
of the one-dimensional reduced action. We show the ex- 
istence of two scaling symmetries and derive the corre- 
sponding Noether currents. In Sec. Ill we discuss the 
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equations of motions and formulate the boundary con- 
ditions of regularity on the polar axis both in the flat 
space and within the full self-gravitating treatment. Us- 
ing them, we derive the sum rules relating the local values 
of some combinations of the metric functions and the YM 
variables to the radial integrals of quantities which ad- 
mit simple physical interpretation in terms of the energy 
density and principal pressures. Sec. IV is devoted to 
study of the asymptotic conditions and investigation of 
the existence of the desired asymptotically vacuum so- 
lutions. Qualitative analysis of the field equations and 
the results of the numerical integration are presented in 
Sec. V. In conclusion we briefly formulate our results and 
sketch further prospects. 



II. THE MODEL 
A. Einstein- Yang-Mills action 

We consider the pure SU(2) Einstein- Yang-Mills action 

1 _ 1 



S = J dx \f—g 
with the field tensor 
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Here e is the gauge coupling constant, and 



A — T A a 



(1) 



(2) 



(3) 



is the matrix- valued gauge potential, with the group gen- 
erators T a normalized as 



T a Tf, = -£ a fc c T c — -^ab- 



(4) 



Variation of the action ([TJ) with respect to the metric 
leads to the Einstein equations 



Ruu — 8vrG ( T av — -g^uTx 



with the stress-energy tensor 



which is traceless: 



pa pat 
fj,a v 



n = 0. 



(5) 



(6) 



(7) 



Variation with respect to the gauge field A® gives the 
Yang-Mills equations 



(D ll ( y /=gFn)* = 0, 



(8) 



where D u is the total covariant derivative with respect 
to both the gauge group and metric connections. 



B. Ansatz for the metric and the gauge held 

We shall use the cylindrical coordinates x^ = 
(t, r, tp, z). The space-time generated by the static cylin- 
drically symmetric source is diagonal and possesses three 
commuting Killing vectors dt,d v ,d z . Its explicit form de- 
pends on the chosen gauge, here we will use the Kasner 
gauge g rr = -l: 

ds 2 = N 2 (r)dt 2 - dr 2 - L 2 {r)d(p 2 - K 2 (r)dz 2 . (9) 

The corresponding non-zero components of the Ricci ten- 
sor are: 



Rn = 



(N'LK)' 
NLK 

(NL'K)' 
NLK 



Dr N" L" K" 



R 



{NLK') 1 
NLK 



(10) 



To construct an appropriate ansatz for the gauge field 
we have to ensure that the effect of the field transfor- 
mation upon the coordinate translation along the Killing 
vector K, 



A a - eC K A 



(11) 



can be compensated by a suitable gauge transformation 


A u — > A u + eD^W, (12) 
with W being a Lie algebra valued function such that 



CkA^ 



(13) 



With three Killing vectors we will thus have for the SU(2) 
YM field nine independent functions. Imposing as usual 
for the static axially-symmetric configurations a discrete 
symmetry M xz ® C, where the first factor stands for re- 
flection in the zz-plane and the second is charge conjuga- 
tion, we reduce the number of independent components 
to six. There are different gauge equivalent ways to pa- 
rameterize the field potential, we choose the following 
gauge: 



eA t = T r wl + T z w%, 

eA r — TipW-^ , 

eA v = T r w^ + T Z W2, 

eA z = T^uig, 



(14) 
(15) 
(16) 
(17) 



where the rotated generators T r ,T v and T z are used: 



T - 1 — 
2,1 



e r = (cos vLp, sin vip, 0), 



(— sin vcp, cos uip, 0) , 



(0,0,1), (18) 



r n being the Pauli matrices. The real number v must be 
integer to avoid multivaluedness of the potential. Two 
temporal components in P^HTT)) can be regarded as elec- 
tric potentials, while the other four are the magnetic 
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ones. Thus purely magnetic static axially symmetric con- 
figurations can be described by four real functions of two 
variables wf (r,z): [IE!, [l6| : 

eA^dx^ = T v (w 2 dr + w\dz\ F (T r w}> + T z vu^ dip. 

(19) 

The corresponding gauge fixing condition can be chosen 
in the form 



ret 



d z w\ = 0. 



(20) 



This ansatz has a residual form-invariance with respect 
to the transformation 



U = exp (T„A(r, z)) 



A 



cos ■ 



A 

2T tp sm-, 



(21) 



under which the functions (u^,w|) transform as 2- 
dimensional gauge field 



d a X, 



a = 1,3, 



(22) 



while the functions {w\ , u>| — v ) behave as scalar doublet: 



cos A sin A 
— sin A cos A 



(23) 



Now we impose the translational symmetry along the 
polar axis assuming that all gauge functions depend only 
on radial variable: wf = wf(r). Performing the gauge 
transformation 



C/ = cxp (T ¥> A(r)) 



tan A 



(24) 



we can exclude the component w\. Furthermore, the 
radial component of the gauge field vanishes by virtue of 
the YM equation wf + X' = 0. Therefore in the static 
cylindrically symmetric magnetic case we are left with 
only two independent functions: 

eA^dx* = r T 2 wjdx 3 + T 3 wldx 2 . (25) 

It is worth noting that with this parametrization the 
equations of motion are invariant under discrete trans- 
formation interchanging both the space and color indexes 
ip <-> z. This symmetry, however, is broken by the bound- 
ary conditions on the polar axis which is different for w 2 
and w\. Renaming the functions wf(r) as in the cosmic 
string theory, we rewrite the final ansatz in the form 



eA^dx 1 " = T v R{r)dz + T z (P(r) 



is) dip. 



(26) 



For another way to arrive at this representation see [17| . 

The corresponding matrix-valued field tensor has the 
following non-zero components: 



T V R' , 



eF in 



T r PR. (27) 



Subsituting them into the action ([T]) we obtain 
S = 



d 4 xy 
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p/2 



It 2 



P 2 R 2 \] 
L 2 K 2 )\ 
(28) 



The stress-energy tensor ([6]) for the gauge field is diagonal 
and can be expressed in terms of the energy density and 
the principal pressures as follows: 



Tq P ^ip F £ 2 "F £ M>) 



rpr 

1 r 

rptp 

<P 



-Pr : 
'Pip 



^ip t-z *W i 



T z z = - Pz = -e 



where 



R r - 



£ip 2e 2 K 2 ' ez ~ 2e 2 L 2 ' ^ " 
Its tracelessness implies 

P=Pr+P v +Pz- 



P 2 R 2 
2e 2 L 2 K 2 ' 



(29) 
(30) 
(31) 
(32) 



(33) 



(34) 



C. Symmetries of the reduced action 

It is convenient to pass to the dimcnsionlcss quantities 
introducing the scale factor I — e _1 v 'AirG and perform- 
ing the rescalings 



r 
N 



N, 



R-^IR, P^ P, 



l — > r x L, 



K —> K. 



(35) 



Integrating the action over t, z, ip we obtain the one- 
dimensional action 



S = ™&&*i s Xt S X = I Ldr 



2/2 



e 2 l 



(36) 



where Vol(i, z, ip) is the normalization volume, and the 
one-dimensional field lagrangian is 



L = 



NL 

p/2 



N'K' 


L'K' 


NK 


+ LK 


R' 2 


P 2 R 2 


K 2 


L 2 K 2 



(37) 



As was argued above, this lagrangian is invariant under 
a discrete tp «-> z transformation which now reads as 



P 



LI \K 



R 



(38) 



In terms of the rescaled quantities the stress-energy ten- 
sor will still be given by the Eq. (|32p with 



R 



2K 2 ' 



pl2 

2L 2 ' 



P 2 R 2 
2L 2 K 2 ' 



(39) 



The dimensional reduction over t, ip and z leads to two 
continuous symmetries. To make this more transparent, 
it is convenient to use an exponential parametrization for 
the metric and gauge functions: 



N 
R 



e , 
e a , 



L = e 

P = e f: 



K = 



(40) 



4 



In terms of these new variables the lagrangian (|37|) 

(y'X' + v't! + X'C - e 2{ P-^(3 12 - 

e («-0+a(/J-A)^ (41) 



L = -e*+ A+ « 



2(a-£) /2 



contains only three independent linear combinations of 
metric functions in the exponential. Choosing them as 
the new variables: 



(42) 



one finds that a, (3 are cyclic. Therefore the lagrangian 
is invariant under the infinitesimal translations of a and 
(3 which leave invariant r, a and ip: 




only two functions. From the Eq. (|32[) it is clear that 
this corresponds to e z — 0, i.e. P = v (to ensure the 
regularity of the axis ) , in which case we find the Melvin 
solution Q, which we will call here "z-string". In the 
general case, in view of the z — (p symmetry mentioned 
above, we deal with the non-linear superposition of the 
z-string and the </?-string. These are interacting via the 
non-linear terms in the Yang- Mills equations. 

It is easy to see that a single ^-string is singular on the 
axis r = 0. Moreover, within the Einstein-Maxwell the- 
ory a configuration which contains both strings inherits 
the (/3-string singularity. However the self-interaction in 
the non-Abelian theory makes possible the superposition 
of both strings regular on the polar axis. As we will see, 
in such configurations the contribution of the ^-string is 
suppressed near the polar axis, this ensures the regularity 
of the full solution. 



' 5a = 





' 5a 


= £2 


sp = 


£l 


sp 


= 


< dtp = 




5tp 


= -£2 


^ = 







= £2 


, sx = 






= 



(43) 



The corresponding Noether charges are 
1 



Qi 



X' - 2e 2 (^ A )/3' 2 ) , 



or, in the previous terms 

Qi = \nlk 



N' V PP' 



1 



-NLK 



~N 



— o — 



(44) 
(45) 

(46) 
(47) 



These are conserved on shell together with the total en- 
ergy (hamiltonian) , which for the static system differs 
from the one-dimensional lagrangian by the sign of the 
potential term: 



H 





( N'L' 


N'K' 


L'K' 




[nl 4 


NK 


+ LK 




p/2 


R' 2 


P 2 R 2 






"^^ 


L?K 2 



= 0. (48) 



Let us discuss some cases studied previously. In some 
gauge models in curved space-time with cylindrical sym- 
metry, such as cosmic strings, the stress-energy tensor 
satisfies the condition 



rjiO rpZ 

1 — 1 z 1 



(49) 



implying the full Poincare symmetry of the (t, z)-plane. 
In this case gu = g Z z and the metric is described by 



III. EQUATIONS OF MOTION AND 
CONDITIONS ON THE POLAR AXIS 



In view of the tracelessness of the stress-energy tensor 
of the gauge field we have three independent Einstein 
equations which read ([5]) 



(N'LK)' 
NLK 

(NL'K)' 
NLK 

{NLK')' 
NLK 



R 



12 



P 2 R 2 



p/2 

U + lO + L 2 K 2 ' 

p^2 R ,2 p2 R 2 
~U + L 2 K 2 ' 

P^_R^_ P 2 R 2 

T? "k 2 ~ Uk 2 ' 



(50) 
(51) 
(52) 





NL\ 










K ) 






NK\ 


N 2 


(, 




= PR 2 




L J 


LK 



In the Kasner gauge © the (rr) Einstein equation is a 
first order equation which is the equivalent to the hamil- 
tonian constraint (|48|) . 

The Yang-Mills equations have a symmetric form 



(53) 
(54) 



As a whole we have the system of five non-linear second- 
order differential equations (I50l - l54|) . which must be com- 
pleted by the set of boundary conditions. In this paper, 
we will be interested in the globally regular asymptot- 
ically Ricci-flat solutions. So we have to formulate the 
regularity conditions on the symmetry axis and to ana- 
lyze the asymptotic behavior compatible with the ansatze 
for the metric and the YM field. The regularity condi- 
tions on the polar axis will be discussed separately for 
the flat-space problem and within the full self-gravitating 
treatment. 
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A. Expansion near the polar axis in flat space 

In the flat space-time N = K = 1, L = r, and the 
equations of motion (|53")l , ([54")) read 



r(rR')' = RP\ 



(55) 
(56) 



For the system of two second-order differential equations 
we have to impose four boundary conditions on the polar 
axis. Requiring the solutions to be regular on the polar 
axis, we find the conditions PR = at r = 0: 



P(0) = fc, R(0) = 0, 



(57) 



where A: is a real parameter. We will consider fc =/= 0, 
since for k = and i?(0) = const the solution is trivial 
R(0) = P(0) = 0. 

Note that in the Abelian case the right hand side of 
the equation ((53)) is zero. Integrating this equation one 
finds a singularity R'(r) = R'(0)/r which can be avoided 
only by excluding the ip— component of the gauge field 
from the ansatz. In the non- Abelian case this singularity 
is avoided by imposing the above conditions. 

To proceed further, we look for the power series solu- 
tions for R and P near the polar axis. The first coefficient 
in the right hand side of the equation for R is propor- 
tional to k 2 . Thus in the left hand side of the equation 
the leading term must be proportional to r k . One can 
find that the other terms in the expansions for R and P 
are proportional to r 2 fe + m ; where n and m are integers. 
In the lowest order one obtains the following expansion: 



R 



P 



qpk 



„2+k 



[fc]! 4(fc + l) 



0(r 



2+2k\ 



pr 



r 2+2k (3( 7 ,3+2fc 



4(fc+ 1) 



(58) 



(59) 



where [fc] denotes the integer part of fc. The free pa- 
rameters q and p represent the boundary conditions to 
the derivatives on the axis r = 0. Together with ([57)) 
one has four boundary conditions which are sufficient to 
formulate the boundary problem. 

For an integer k the parameters q and p are equal to 
the higher order derivatives of R and P as follows 



p = P"(0), 



R {k) {0), 



(60) 



and the expansions contain only integer powers of r, 
therefore the solutions for R and P are analytic func- 
tions. In this case we can obtain the higher order terms 
in series of expansions which are presented in Tab. |I] 
These have the following structure: 



r k+2j 



1=0 



p = k+ P r -+Y^p, 



2(fc+j + l) 



3=0 



r 2(k+j+l) 

W+TTTy. 



(61) 



, (62) 



where the coefficients of any order can be expressed 
through the free parameters p, q as follows 



R. 



k+2] 



[HI J 



i=0 



(63) 



P 



L k + 1 J 



2(k+j+l) 



i=0 



#g2(.'+V'-*+l). (64 ) 



Thus we have the three-parameter family of solutions 
to the system (|551l56|) specified by the set {q, p, fc}. But 
as will be shown below one can fix e.g. q = 1 and con- 
sider only the two-parameter family {p, fc}, while other 
solutions with q ^ 1 can be obtained by transformations 
of the family 5=1. 

It is easy to see that if {R(r), P(r)} is a so- 
lution to the system ([551 EH) then {lR(lr), P(lr)} 
will be also a solution, where I is some con- 
stant. The question arises whether this rescaling 
affect parameters {q, p, fc}. Consider the solu- 
tion {R(qo,po,ko;r), P(qo,po,ko;r)} and make the 
scale transformation {lR(qQ,po,ko]lr), P(qo,po, fco;Zr)}. 
The latter solution must coincide with some solution 
{R(qi,Pl,kf,r), P(qi,Puh;r)}. In the leading orders of 
expansion one has 



qir k = lq r k l k 7 
h+ Pl r 2 /2 = k Q +p r 2 l 2 /2. 



(65) 
(66) 



This implies that fc; = fco = fc, pi — pi 2 and I = 
{<li / 1o) TTT ■ Extracting I one obtains 



P = Po(<l/qo) k+1 - 



(67) 



Thus all solutions belonging to the curve (|6T|) in the plane 
fc = const of the parameter space can be obtained from 
the solution {go = 1, Po, fc} only by the scale transfor- 
mation 

R(q,p,k;r) = q^TT R (l,pq~^TT, fc; q^rj , (68) 
P(q,p,k;r) = P (l,pq~*^ , fc; q^r^j . (69) 
In the following we will restrict to the case q = 1- 

B. Curved space 

Near the polar axis the metric component g w = L 2 (r) 
must vanish. To ensure the finitcncss of the curvature the 
regularity conditions on the metric functions and their 
derivatives must be imposed. To find them we consider 
the Kretchmann scalar, 



N'L' 
NL 

N 



N'K' 



L 



NK 

2 



TF I + -T + hr , (70) 



K" 
K 



L'K' 
LK 

2 
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which must be finite on the polar axis. Since L vanishes 
there, one has to require 2/'(0) — as well. This implies 
that 2/(0) must be non-zero, otherwise g vv will be iden- 
tically zero. Furthermore, vanishing of denominators in 
the upper string of (|70|) must be compensated by van- 
ishing of numerators, which implies N'(0) = K'(Q) = 0. 
Finally, N and K cannot vanish on the axis or they will 
be equal to zero everywhere. Therefore the locally flat 
metric starting with 



Comparing this expression with the component F rtp of 
the field tensor (|27|) . one can see that this quantity is 
nothing else than the corresponding tetrad component 
with the color index z. It is therefore the limiting value 
of the physical magnetic field component along the sym- 
metry axis 



B z (r) = - 



P'(r) 
eL(r)' 



(76) 



N(0) = K(0) = 1, 
AT'(O) = K'(0) = 0, 



L(0) = 0, 
2/(0) = 1. 



(71) 



is the only one (up to rescaling) which provide finiteness 
of the expression ([70)) one the polar axis. Note that from 
the azimuthal Einstein equation one can show that the 
condition 



After rescaling, we obtain the dimensionless parameter 

1 = 47rGe" 1 B(0), (77) 

therefore the parameter p is equal up to sign to the (di- 
mensionless) value of z-component of the magnetic field 
on the polar axis: 



L"(r = 0) = 



(72) 



is satisfied if the azimuthal pressure is finite near the 
polar axis. 

Now we use the same strategy to find the appropriate 
condition on the gauge field (|2"6")l constructing the power 
series expansions near r = for all variables satisfying 
the metric conditions (|7T|) directly from the system (J5DJ- 
[54l) . One can expect similar results for the field functions 
as in the flat case with some amendments due to gravity. 
This is indeed the case as can be seen from Tab. [TTJ But 
the expansions of the metric functions for non-integer k 
read 



P = 



Sum rules 



(78) 



To relate the expansion on the symmetry axis to the 
asymptotic form of the solution on can derive the sum 
rules following the approach of the Ref. [l8[ . Multiplying 
the Einstein equations ([HO] IE2]) by 



-9 



NLK, 



(79) 



N ~ l + C N r 2k + 0{r 2k+2 ), 
K ~ l + C K r 2k + 0(r 2k+2 ). 



(73) 
(74) 



we integrate them from r — to r using the expressions 
for the stress-energy tensor. This gives 



Obviously the terms in the Kretchmann scalar (|7D|) which 
are proportional to the second derivatives of N and K 
vanish in vicinity of the axis only when k > 1. 



„2fc-2 



N'LK = M, 
NL'K = W+l, 
NLK' = X, 



(80) 
(81) 
(82) 



The energy density and pressure, which are proportional 

to r 2k ~ l in the vicinity of the axis, vanish in this case. where three new functions are introduced to the integrals 



Therefore the condition ([72]) is also satisfied. 

Note that the gravity coupled solution does not contain 
any new free parameters. The full set of ten boundary 
conditions consists of the previously obtained two condi- 
tions for the field functions ([57]) and the six conditions on 
the metric (|7ip . with three parameters {p, q, k} remain- 
ing free (recall that to ensure analyticity of the solution 
near the polar axis one has to take integer k). 

Unlike the flat case, now the series expansions contain 
some additional terms. But all such terms include only 
higher powers of parameters p and q apart from the main 
terms from (|58p. (159")) . For this reason in the weak field 
case p — > 0, <7 — > the curved space expansions are dom- 
inated by the flat space ones ([35} , ([59"|) , i.e. the solutions 
satisfying ([fJ7|) will differ only by the scale factors. This 
scale invariance does not hold in general, however. 

Now we would like to clarify the physical meaning of 
the parameter p. Using the Tab. HTlone finds 



M = 



W 



X 



dr V-g(P + Pr + P V + Pz), 


r 

- J dr</=g(p - p r + p v -Pz), 
drV^g(p ~ p r - Pip +Pz)- 



(83) 
(84) 
(85) 



P 



lim 



P'(r) 
6 L(r) ' 



(75) 



The function M(r) is the Tolman mass of the gravitat- 
ing YM field within the cylinder of the radius r. Two 
other quantities are related to integral pressures in the 
longitudinal and azimuthal directions, their more precise 
meaning will be clear later on. The sum rules (|80H82j) 
express the metric functions through the integral charac- 
teristics of the matter field. 

Two other some rules can be obtained by integrating 
the YM equations. Multiplying (J33J [5j} on R and P 
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respectively, we can rewrite the YM equations as follows 



/ NL 
NK „„, 



R' 2 


P 2 R 2 


K 2 ' 


h L 2 K 2 


p/2 


P 2 R 2 


U ' 


h L 2 K 2 



XLK ! —7 - -T-r— 7 i ■ (86) 

(87) 



Then we integrate in the same way as the Einstein equa- 
tions obtaining the sum rules: 



NL (R 2 
~K 2 



_ . R' 2 P 2 R 2 



NK {P 
~L 



2y 



P' 

lim I NKP— 

r^O \ L 



M-X 



P 2 R 2 



pK 

L 2 ' L 2 K 2 
M-W 



m . (89) 



From the three integral functions M, W, X only two are 
independent: the constraint equation (|48[) implying 



(90) 



MX + M(W + 1) + X(W + 1) = 0. 



IV. ASYMPTOTIC BEHAVIOR 

Methods of investigations of the asymptotical geome- 
try generated by cylindrically symmetric field configura- 
tions were developed within the model of Abelian cos- 
mic strings @, [H, H3, [Sj], f° r a more general model- 
independent approach see [HI]. Here we shall use the 
same technique with some modifications appropriate to 
the YM field dynamics. 



Kasner solution 



integral functions M, W, X assuming them to be finite 
Moo = lim M(r), 

r — >oo 

Woo = lim W(r), (93) 

r — >oo 

Xoo = lim X(r), 



and denote £oo = Moo + Woo + X c 
parameters satisfy a constraint (f9"0"l) : 



1. These three 



MooYoo + Moo(Woo + 1) + X^W^ + 1) = 0. (94) 

Together with B, we thus have four parameters to be 
related to the Kasner parameters. From the Eqs. Q80J - 
1521 it follows 



N'LK 
NL'K 



Mo 



Woo + 1, 

NLK' -» Xoo- 



(95) 
(96) 
(97) 



It is easy to find the solution to this equations assuming 
Soo ± 0: 



N 
L 
K 

-g 



cur a , a = Moo/Soc, (98) 

7^ b , & = (VFoo + (99) 

(^ 7 )"V C , C = Xoo/Xoo, (100) 

YiA- = E^r, (101) 



where cj,7 are some constants. The first Kasner con- 
dition (f92|) is obviously satisfied , while the second is 
equivalent to (|M)l . 

So we are left with only two independent parameters 
describing the asymptotically vacuum metric. The condi- 
tion ([M]) means the fulfilment of the radial Einstein equa- 
tion. From (jlOip it is clear that the condition ^ is 
the requirement of the non-degeneracy of the metric at 
infinity. On the polar axis r = the boundary condition 
for y/— g read that its derivative is positive. Since y/—g 
must be non-vanishing at any point except r = 0, one 
has to require > 0. 



It is well known that the general vacuum solution of 
Einstein equations for the ansatz ([9]) is given by the Kas- 
ner metric 

ds 2 = r 2a dt 2 - dr 2 - ~f 2 r 2b dtp 2 - r 2c dz 2 , (91) 

where constants a, b and c satisfy the conditions 

a + b + c = a 2 + b 2 +c 2 = 1. (92) 

Consider the cylindrically symmetric configurations of 
the YM field for which the stress-energy tensor vanishes 
at r — > oo. Then the metric in the asymptotical region 
will approach (|9"T|) , where parameters a, b and c must be 
specified in terms of parameters of the field configuration. 
For this purpose we define the asymptotic values of the 



B. Existence of asymptotically vacuum solutions 

Now let us turn to the dynamic of the gauge field. 
The assumption of the finiteness of Moo, Woo and Xoo 
imposes some conditions on the gauge field asymptotics. 
Actually, the energy density must decrease fast enough 
for Moo to be finite. From Eqs. ((551 -155 11 it is clear that 
the absolute values of Woo and Xoo are smaller than Moo ■ 
The total energy density is the sum of three components 
((29l) . which asymptotically can be written as 

e v ~i?' 2 r- 2c , e z ~P' 2 r- 2b , e w ~ P 2 R 2 r 2a - 2 . (102) 

All these quantities must fall faster than r~ 2 to ensure 
finiteness of Moo- Notice that 0<a =M 00 /E 00 <1 since 
Soo > 0, and Moo > because it is proportional to the 
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energy of the configuration per unit length. Therefore 
2a - 2 > -2 and RP = o(l) as r — ► oo. 

Now we can obtain further information invoking the 
sum rules ([881 I89[) in combination with the metric asymp- 
totics (P HTUD|) : 



(R 2 y 
(p 2 r 



(Moo - X^)K 2 

(Moo^Woc_- TO ^ 2 



(103) 
(104) 



Integrating these equations in the leading order with ac- 
count for igSHTTJOl) we obtain 



R 2 = C {R) 

P 2 = d p) 



C { 2 R) T 2C , 



(105) 
(106) 



with four constant parameters, which can be determined 
using ([102)1 . The integrand in (l88|) is positively defi- 
nite, so the integral is non-zero. Since C% ^ s propor- 
tional to this integral, it must be non-zero too. To sat- 
isfy the fmiteness condition for with asymptotic be- 
havior of the first exponents in Eqs. (|102p the quantity 
2(2c — 1) — 2c = 2c — 2 has to be smaller than —2, so c 
must be negative. Since among the set of parameters a, b 
and c satisfying the Kasner conditions there can not be 
two negative quantities, b has to be positive and the cor- 
responding exponent in Eqs. ([102|) is equal to 26— 2 > —2. 

Thus to ensure the finiteness condition for M m we must 

(p) 

take Co =0, or, in other terms, 



p/2 p2 7->2 > 

drNLK | — + -777-7777 I - m = 0, (107) 



L 2 L 2 K 2 



what is equivalent to one more constraint on the param- 
eters: 



2kB. 



(108) 



Note that a similar relation (up to some factor) was ob- 
tained within the model of Abelian cosmic strings [2l[ as 
the condition of existence of the Melvin-like solutions. 

To get further information about asymptotics of the 
gauge functions we apply the same integration procedure 
as described in section III to the non-transformed equa- 
tion of motion ([53"[) . This gives 



R = d+ C 2 r 2 



(109) 



with two more constants. Since c < 0, for the square 
roots of (|105|) one has two possibilities: 



Comparing with (|109|) we see that the second case is re- 

(P) 

alized. But this implies that C\ =0, and P — > 
at infinity faster than power-law, while R approaches a 
constant as a power with the exponent 2c. Finally one 



finds the following explicit expression for C, 



OR). 



a 



(R) 



4C 1 X 00 w 2 7 2 • 



(112) 



Now consider the flat space case. We are looking 
for asymptotically vacuum configurations, which implies 
vanishing energy density at infinity: 



rR'- 



p/2 

r 



R 2 P 2 



0. 



(113) 

Integrating the equations of motion (j55[) , (156|) with zero 
right hand sides, one obtains the asymptotical expres- 
sions for field functions 



R = Cilnr + Ci, 
P = C 3 r 2 + C 4 , 



(114) 
(115) 



which satisfy ([113[) if C3 = 0. Also there must be C4 = 
otherwise the right hand side in the equation for P will be 
non-zero. Thus the conditions for the field components 
at infinity read 



R~lnr, P^0. 



(116) 



These are the necessary conditions of vanishing of the 
energy density. For more detailed expansion see sec.V. 

Performing the same steps with the equation for P, 
one obtains the following relation: 



(P 



2y 



lim 

r — >oo 2r 



drr 



P' [ 



R 2 P 2 



P' 

lim P— . (117) 

1 — >o r 



Therefore the condition (| 1 16[) is satisfied if the integral 
exists and the whole expression is equal to zero. This is 
the necessary condition of the existence of asymptotically 
vacuum solutions: 



drr 



pi. 



R 2 P< 



(118) 



In the flat case all values M(r), W(r) and X(r) diverges 



as lnr at infinity because of the term R jr 



the 



integrand. But the difference (M — W) remains finite 
and it is equal to the integral in the left hand side of 
Eq. ([118[) . So the form of the bifurcation condition (1108)) 
remains the same though the separate values of and 
Woo are not well defined. 



R = 
R = 



c[ R) r c , 



a 



(R) 



(R) 



0, (110) 

r 2c , C[ R) f 0. (Ill) 



C. Kasner parameters 

So far we have obtained the asymptotic behavior of the 
gauge and metric functions consistent with the equations 
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of motion. The solution for the metric is determined by 
three parameters a, b and c satisfying two constraints. 
Now recall that there are also two conserved Noether 
charges (|46p , P?} . This enables us to expose the asymp- 
totic parameters of the solution through the data on the 
polar axis. 

Equating the Noether charges Qi\ r =o = Qi\ r ^oo onc 
obtains further two constraints: 



1 = £oo(a-6), 



= Soo a - c - 



Moo — 



(119) 
(120) 



These relations can be used to express b and c in terms 
of two new parameters: 



b = a + P, (3 = 



a — cr, a = 



1 - 2fcB 

Moo ~ 



(121) 
(122) 



Now, from the first Kasner condition one has 



The parameter v does not enter the equations of mo- 
tion. This parameter must be integer to avoid the mul- 
tivaluedness of the solution. Actually it (and only it) 
specifies the azimuthal part of the ansatz, while k (and 
only k) specifies the radial part. The full analytic solution 
comes as a product of both parts and thus is character- 
ized by two integers. The natural question arises whether 
the integers v and k can be given any topological inter- 
pretation. 

Using the asymptotic expansion for the gauge potential 
A M (recall that for large r the function P vanishes while 
R tends to a constant, C), we find asymptotically the 
potential one-form as 



A\ oa = CT v dz ~ vT z dip, 



(127) 



and the field two- form F = dA + A A A vanishes. This 
implies that A in (| 1 2T[) describes a vacuum state and can 
be expressed through the matrix-valued gauge function: 

A\ x = UdU~ l . (128) 

This equation can be solved as follows: 



a=-(l + a-0). 



(123) 



Finally, the quadratic Kasner condition allows us to ex- 
press all parameters in terms of the only parameter /3: 

1 



a = -(1-309/2) + VI -309/2)'), (124) 



b = -(1 + 309/2) + VI -309/2)8), (125) 



c = -(1-2V1- 309/2)2). 



(126) 



In the previous subsection we have showed that c must 
be negative, this determines the sign of square root, an 
alternative choice of a sign would give c > 1/3. We will 
show later that the expression under the square root is 
always positive. 

Thus the value B of the magnetic field value on the 
symmetry axis along with £m/t y fully specify the solu- 
tion in the asymptotic region. 



D. Asymptotic YM vacua 

As was shown above the real number k specifies the 
family of solutions for R and P i.e. the radial part of 
the ansatz. Since R ~ r k in the vicinity of the axis, the 
azimuthal component of the magnetic field B v ~ R' ~ 
r fc_1 diverges for non- integer k on the interval < k < 
1 and so does the Kretchmann scalar ([70]) . To ensure 
analyticity of the solution near the polar axis k should 
be chosen integer. Thus k can be regarded as a "radial" 
number of the configuration. If k = there are no regular 
solutions except the trivial one, since the condition (|108|) 
implies vanishing of the energy density. 



f7=-^=(cos/(l-2T r )-2sin/(T ¥3 + T 2 )), (129) 



where 



f=(Cz- u<p)/2. 



(130) 



Now consider the potential on the polar axis. Here it 
is also a pure gauge: 



Al^s = (k- u)T z d<p = UdU- 1 



(131) 



with 



U = cos/- 2sin/T z , f=(k-v)<p/2. (132) 

These two vacua are related by the gauge transformation 
(I129P with / = (Cz—k(p)/2. The winding number of this 
transformation is zero and the transformation is "small" . 
This reflects the fact that the mapping of the asymptotic 
cylinder (r — * oo, z) into the gauge group SU(2)^ S 3 is 
always topologically trivial. Thus we can conclude that 
the integers v and k do not have any topological mean- 
ing. Note also that the divergence of the Chern-Simons 
current 



lbn' s 



trFa V F^ 



(133) 



is zero because of the absence of electric part in the 
ansatz. 

It is instructive to compare the situation with that in 
the case of the EYM sphalerons of Bartnik-McKinnon's 
type. Sphalerons are the finite energy saddle points in 
the configuration space of the system, whose existence 
can be revealed by the constructing the non-contractible 
loops of configurations [22j which alternatively can be 
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seen as connecting topologically distinct vacua. For the 
spherically symmetric four-dimensional EYM configura- 
tions this was done explicitly in [23]. Configurations of 
finite energy along the loop form a manifold homeomor- 
phic to S 3 , and the existence of non-contractible loops 
is a consequence of non-triviality of the mapping of this 
S 3 to the gauge group SU(2) ~ S 3 . In the cylindrical 
case one considers the loop of configurations of finite en- 
ergy per unit length which form the manifold equivalent 
to S 2 . Thus instead of the mapping S 3 — ► S 3 one deals 
with the mapping S 2 — > S 3 which is topologically triv- 
ial. Therefore in the cylindrical case the solutions do not 
admit sphaleronic interpretation. 



V. EXISTENCE OF GLOBAL SOLUTIONS AND 
NUMERICAL INTEGRATION 

In the previous sections we have shown that global 
solutions asymptotically approaching vacuum configura- 
tions, if exist, should be fully specified by the boundary 
conditions on the polar axis. Up to rescaling, the only 
relevant quantity is the parameter p = —M which is the 
second derivative of P on the axis (see Sec. Ill, we set 
the scale taking q = 1). It turns out that the family 
of solutions specified by different p exhibits the bifurca- 
tion structure similar to that observed in the Bartnik- 
McKinnon case, so the existence of regular solutions can 
be proved by similar methods. Here we sketch a proof 
and present the results of numerical integration which 
strongly supports this analysis. 



A. Flat case 

To better understand the dynamics of the gauge field 
we start with the flat-space problem. Our goal is to in- 
vestigate the bifurcating behavior of the system ([55l l56|) . 
Unfortunately the dynamical systems technique does not 
provide enough information about solutions. There are 
three special points of the system: two at the the origin 
(r = 0, P = 0) and (r = 0, R = 0) and one at infinity 
(1/r = 0, P = 0). The corresponding linearized systems 
are degenerate and phase portraits in the vicinity of the 
special points are trivial. It turns out that the values 
of parameters p and q are non-small on the bifurcation 
curve, so we can not use the asymptotic methods for the 
analysis. Our reasoning will based on the fact that the 
equations (|55l [56)) are linear with respect to the variable 
entering with its derivative, if the other variable is con- 
sidered as a given function. 

First, let us investigate properties of the solutions of 
the Eq. (|55[) with an arbitrary regular function P(r). In 
what follows we assume P(0) = k > 0. 

Lemma 1. The solution for R(r) with the boundary con- 
ditions R(0) = 0, R'(0) = q > is monotonously increas- 
ing function for r > 0. 



Proof. Since R'(0) > 0, it remains positive in some region 
< r < a. Assuming that it vanishes at r — a, we 
integrate the Eq. (|55|) on the segment [0, a]: 

a 

1 f RP 2 

R {a) = ~ / — dr. (134) 
a J r 

o 

Since R 1 is positive on [0, a), the function R itself is also 
positive, R(r) > for r e (Q, o]. Due to the assumption 
P(0) > the function P 2 (r) is not identically zero on this 
segment. Therefore the integral at the right hand side of 
the Eq. (|134[) is strictly positive and hence R'(a) ^ 0, 
what contradicts to the initial assumption. □ 

Lemma 2. Consider two positive functions satisfying 
P\{r) > -P2M for r > 0. Then the difference AR = 
R[q, Pi] — R[q, Pi] is monotonously increasing function 
for r > 0. 

Proof. Using the expansions near r — one can see that 
the rate of growth of AR near the origin is proportional 
to the value of P 2 — P 2 at the origin. Therefore there 
exists some small value r — a for which AR(a) > 0, 
AR'(a) > 0. Again assuming that at some point r = b 
the derivative AR vanishes and integrating the equation 
for this quantity (which is essentially the same as for R' 
due to linearity of the Eq. [53)1 we obtain 

b 

AR'(b) = AR'(a)y + I [ AR ( P ?~ P ^ dr. (135) 
b b J r 

a 

Similarly to the previous lemma, we obtain AR'(b) > 
what contradicts to the assumption. Thus AR(r) is 
positive and AR(r) is monotonously increasing for r > 
0. □ 

Now let us turn to the properties of the solution of 
the Eq. ([56]) for P with a given arbitrary regular function 
R(r). 

Lemma 3. If at some point r = a one has P{a) > 0, 
P'{a) > 0, then both P(r) and P'(r) are positive 
monotonously increasing functions for r > a. 

Proof. Assume that there is some point b > a such that 
P'(a < r < b) > 0, P'(b) = 0. Then P will be a non- 
decreasing function on the segment [a, b] implying P > 0. 
Direct integration of the Eq. (|56|) leads to 

b 

P'(b) = P'(a)b/a + bj R 2 Pdf. (136) 

a 

The integrand at the right hand side is non-negative, 
therefore P'(b) > contradicting to the assumption. 
Thus P(r) is a monotonously increasing function and 
from the Eq. (ll36p it follows that P' is monotonously in- 
creasing for r > a too. □ 
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Since the Eg. (|56|) is linear in P, it is invariant under 
inversion P — ► —P, therefore if P(a) < 0, P'(a) < then 
— P(r) and — P'(r) are positive monotonously increasing 
functions for r > a. 

Now let us turn to the full system (|5"B1 EU) ■ One can 
easily see that there exist integral curves corresponding 
to both cases considered above. Indeed, for small values 
of \p\ and large values of q the derivative P' with the 
growing of r soon becomes positive with P also remaining 
positive. Vice versa, for large \p\ and small q the function 
P soon changes the sign while P' also remains negative. 
In both cases we will observe an infinite growth of \P\ 
as r — > oo. Between these two sets of p, q there is a 
bifurcation set such that the corresponding integral curve 
satisfies the inequalities P(r) > 0, P'(r) < as r — ► oo. 
Then \P\ will be a monotonously decreasing function at 
infinity. Thus we will obtain the heteroclinic trajectory 
connecting two special points of the system: the origin 
(r = 0,P = 0) and infinity (1/r = 0, P = 0). Using the 
above lemmas one can prove that this solution is unique: 

Theorem 1. For each set of parameters q,k there is a 
unique value of parameter p corresponding to the hete- 
roclinic trajectory connecting points (r = 0, R — 0) and 
(l/r = 0,P = 0). 

Proof. Let us assume that there are two such trajectories. 
Since the solution of the system is fully specified by the 
set of parameters (q,p, k), different trajectories must cor- 
respond to different values of p, say p\ > Pi- The series 
expansion for P near the origin implies AP = Pi — Pi « 
(pi — pi)r 2 /2 > 0. This means that there is some point 
r = a, where AP(a) > 0, AP'(a) > 0. Now consider the 
point r = b where AP' vanishes: AP'(a < r < b) > 0, 
AP'(6) = 0. Since AP' is positive, the function AP is 
monotonously increasing which implies that Pi > Pi on 
the segment [a, b] ( with both Pi, Pi remaining positive). 
Using the lemma 2, we find that Pi > Ri on this seg- 
ment. Integration of the equation for AP gives 

b 

AP'(b) = AP'{a)b/a + b J (P?P X - R\P 2 )df. (137) 

a 

The integrand is positive and the quantity AP'(a) is pos- 
itive too, therefore AP'(6) > what contradicts to the 
assumption. This means that AP(r) is monotonously 
growing function for r > a with the rate of growth not 
slower than r 2 . But as r — ► oo, both Pi, Pi vanish and 
AP has to vanish too. Thus we arrive at the contradic- 
tion again, hence the assumption about the existence of 
several heteroclinic trajectories is wrong. Therefore there 
is a unique value p = pb for each set (g, k) for which the 
corresponding solution satisfies the boundary conditions 
P(0) = 0, P(0) = k 7^ and monotonously reaches its 
asymptotic form (| 1 16|) . □ 

Thus the full solution space consists of the following 
sets: 



• Solutions with monotonously increasing |P|. Nu- 
merical calculations show that these grow infinitely 
at some finite r, i.e. they are singular. 

• Solutions with monotonously decreasing | P|. There 
is a unique bifurcation value p = pb for each set of 
(q, k) for which such a solution can arise. These 
solutions are globally regular. 

Numerical integration of radial equations starting from 
the polar axis r = shows that for any fixed value of q 
both types of solutions are realized. Singular solutions 
diverge at finite distance r = r s as follows: 

R(r) -> +oo, P(r) -> +oo, as r -> r s , (138) 
R(r) — > +oo, P(r) — > — oo, as r — > r s . (139) 

In the vicinity of the singularity, the system (|55| , ([56)) 
can be integrated analytically. In the leading order one 
has 

R" = RP 2 , (140) 
P" = P 2 P, (141) 

where P = P/r s . This system is invariant under the 
transformation R <-> P (actually this is the symmetry 
(|38f V Multiplying these equations on P and P respec- 
tively one obtains 

PP" = PP," (142) 

and thus there exist solutions P = ±P. Then the equa- 
tion 

R" = P 3 (143) 

can be easily solved, and find we the divergent terms as 
follows: 

P=^-, P = ±^i. (144) 
r s - r r s - r 

Now discuss globally regular solutions. It was con- 
firmed numerically that for any fixed value of q there is a 
unique value p = Pb such that for any p > pb and p < pb 
the integral curves correspond to singular solutions of 
the first and second types respectively. As p — > pb the 
singularity point r s monotonously moves away from the 
origin presumably to infinity. In order to reconfirm the 
existence of regular solutions numerically we have per- 
formed integration from infinity gluing it with the so- 
lution started from the polar axis. Asymptotically the 
regular solution to the system (|55H56| is given by the 
following three-parameter family: 

R= (Cilnr + C 2 )a{r), (145) 

P = dJ T ^exp(-C 1 r]nr+(C 1 -C 2 )rW(r), (146) 
V mr 

where some functions a(r), (3{r) satisfying the boundary 
conditions a, (3 — >lasr— >oo. The expansions for 
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a, (3 are rather complicated because of the presence of 
logarithmic terms R and exponential terms for P. In 
particular, 



x = r 



-1/3. 



0{r) 



C 2 



2Ci lnr 



+ o(ln r), 



(147) 



a(r) = l+(^ 7 ) 2 (l + «(ln- 1 ? -))- (148) 



Therefore any globally regular solution is fully speci- 
fied by the set parameters (C\,C<z,D) as well as by set 
(q,p,k). In the asymptotic set the parameter D defines 
a scale, reflecting the situation with the set (q,p, k). Nu- 
merical integration provides matching of these two sets 
of parameters with each other. 

Previously we have obtained the necessary condition 
(|118| of existence of regular asymptotically vacuous so- 
lution. Numerical integration confirms that for the so- 
lutions at the bifurcation point p — pt this condition is 
satisfied indeed. Both regular and singular numerical so- 
lutions for q = 1, k = 1 are shown in Fig. [TJ For other 
values of k the behavior of solutions is similar. Also, as we 
have argued there are families of the solutions related by 
the scale transformation. Numerical computations con- 
firms that all solutions with 



P = Pbq l 



(149) 



are regular. They correspond to the bifurcation curves 
of in the parameter space shown in Fig. [2] 



B. Self-gravitating case 

As was argued before, the behavior of the full self- 
gravitating solutions corresponding to small values of pa- 
rameters q — » 0, p — > must be the same as in the flat 
space. This implies the existence of both regular and sin- 
gular solutions in the gravitating case too. Numerical in- 
tegration of the system (|53ti50| proceeds along the same 
lines with the boundary conditions obtained from the ex- 
pansions on the polar axis with free parameters (q, p, k) 
(Tab. II). The main difference is that now the asymp- 
totic behavior of the metric functions and R at infinity 
is power-low with the linear combinations of the Kasncr 
parameters (a, b,c) in the exponent. For different sets of 
these parameters the linear combinations with the same 
coefficients can be positive or negative which mean either 
infinite growth or vanishing of the corresponding terms in 
the expansion. For this reason the asymptotic solutions 
will be different for different sets of (a, b, c). We give here 
the asymptotic expansions in the most symmetric case of 
the values of parameters: a — 2/3, b — 2/3, c = —1/3. In 
this case the functions (N, L, K, R) expand in powers of 



R = d + C 2 x 2 



N 



24Cf x 4 
w 2 7 2 4! 

?3 



+ 0(x 5 ), 



(150) 



x 2 W7 2 7 w 3 7 4 3! 
7 



L = N— 

LU ' 

x 4 

K = lujx-UujCs— + 0(x 5 ). 
The YM function P has an exponential form: 



(152) 
(153) 



P — exp < — 



4to>7 



2C 2 
x 2 



4C * Inx) \{D+0{x)). 



(154) 



w 2 7 2 



Altogether we have six parameters: Ci, C2 and D for 
the YM functions as in the flat space, and four metric 
parameters: u), 7, C3, two of which can be changed by 
rescalings of i and z, and one independent Kasner param- 
eter. Together with three parameters on the axis g, p, k 
one has ten parameters to match the solutions obtained 
by the integration from the axis and from the infinity. 
Since these solutions come from the system of five second 
order differential equations the number of parameters co- 
incide with the number of matching conditions. The nu- 
merical integration scheme is the same as in the flat case, 
and the conservation of charges Qi, Q2, H along the in- 
tegral curves is used to control the numerical error. 

It turns out that in the gravitating case the bifurcation 
parameters do not satisfy the condition (|149[) anymore. 
The bifurcation curve has to be found numerically. Ac- 
tual calculations show that it is bounded within the finite 
region 



< q < q n 



0>P>Pn 



(155) 



and for each q in this region there are two bifurcation 
parameters p. The shapes of the bifurcations curves are 
similar for all k, but with growing k the region of param- 
eters becomes very narrow, this is illustrated on Fig. [3] 
for three lowest k. The condition (|107[) is satisfied in- 
deed for the regular solutions, and all other analytical 
results of the Sec. IV are also confirmed by numerical 
computations. 

Using the Eqs. (|124H126]) we can relate the asymptotic 
Kasner parameters to the value of the magnetic field on 
the symmetry axis. The weak field limit B — > corre- 
sponds to the region q — > 0,p — > on the bifurcation 
curve. In this case the metric approaches the Minkowski 
metric. From the asymptotic relations (f9"8"l - |101[) one can 
see that this implies M^, Woo, -— > 0, Sqo — ► 1, and 
the parameter j3 in Eq. (|124N126j) approaches the value 
1 — 5, S — *■ 0. In the first order in 5: 



N 



K 



(156) 



As we move along the bifurcation curve, the field pa- 
rameter B grows up, the exponent of L decreases, while 
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the exponent of N increases. For B = I /(2k) the metric 
acquires an additional boost symmetry: from Eq. Q12ip 
it follows that (3 — so the Kasner parameters are 
a = b = 2/3, c = —1/3. For this particular solution 
the metric possesses the asymptotical boost symmetry 
in the plane (t, ip) since gtt becomes proportional to g vv 
as t — ► oo . 

With B growing further, the Kasner parameters sat- 
isfy a > b. But still gtt and g vv grow at infinity while 
g zz decreases. For the Melvin universe one has growing 
9zz — 9tt an d decreasing When the bifurcation curve 
approaches the p-axis, which corresponds to the Abelian 
case with R = 0, the metric is of the Melvin type. It is 
interesting to observe how these different types of metric 
transform into each other. It turns out that for q — ► 
but finite p the functions R, R' are both small, while 
P, P' are not. This is why only the z-string (associated 
with the function P) persists in the vicinity of the sym- 
metry axis. It generates the Melvin metric with growing 
r, while the energy density decreases as a power-law. In 
this region the stress-energy tensor exhibits the Melvin 
symmetry: 



rpO 

1 o 



rptp 
V 



(157) 



But for larger r, the functions R and R' become non- 
negligible, and the energy density of the tp— string in- 
creases and reaches its maximum Fig. [5] In this re- 
gion the energy density of the 2— string decreases expo- 
nentially. So at the small distances from the axis the 
z— string dominates and generates the Melvin metric, 
while at larger distances the tp-string dominates, trans- 
forming the metric asymptotically into the form ([124I - 
I126p . In this region the components of the stress-energy 
tensor satisfy the relations 



1 o 



TV 



(158) 



As p — > Pmax , the boundary between these two regimes 
moves to infinity. Since strings cores are largely sepa- 
rated, one can introduce the effective radius ro such that 
the mass of the z — string M z is concentrated in the re- 
gion (0, ro), while the mass of the ip— string M v corre- 
sponds to the region (ro, 00). We define then 



M z = M(r ), M v 



M(r ), (159) 



and similarly for the quantities W, X, E. The solution 
within the region (0,ro) can be considered as Melvin, so 
it is known that M z = 2, and due to the symmetry (|157[) 
one has M z — —W z = X z = 2, £ z = 3. The integration 
of the Eqs. JMUZD gives: 



N'LK 
NL'K 
NLK' 



M v 



-M z 

-w z 

■X* 



■■ M v + 5 
1 = W w 



X^ + X z — X v + 2. 



(160) 
(161) 
(162) 



Using the asymptotic symmetry (|158|) . one finds 



in terms of the ^—string mass M v only: 



M v + 2 
M w + 3' 



1 



M,„ + 3 1 



2-M v 
M v +2>' 



(163) 



For them the first Kasner condition is satisfied automat- 
ically, and the second gives the following equation for 



M V (M V -A) = 0. 



(164) 



So we have two possibilities: either M v — implying 
the pure Melvin metrics, or M v = 4 (in our dimension- 
less units) which gives the limiting solution a = 6/7, b = 
3/7, c = —2/7. For small values of the parameter q — > 
when the above picture of well separated strings is valid, 
the mass of the ^—string approaches the value 4. Re- 
call that the pure Melvin case, the mass of the z— string 
approaches the value 2 as p — * 0. 

The numerical results for k = 1 are presented on Fig. [U 
where the dependence of Kasner parameters on the field 
strength p is shown. A typical numerical solution is 
shown on Figs. 15161 The Fig[7] describes how the com- 
ponents of the stress-energy tensor switch from the sym- 
metry (|157|) to (|158|) with growing r. On Fig. [5] one can 
see the core of the z— string in the vicinity of the axis 
which forms the Melvin metrics. Fig. [5] demonstrates 
that at some distance from the axis there is another pick 
of the energy density — the core of the ip— string, which 
changes the signs in the exponents of g zz and g vv . 

Since at p = p max both Abelian and non- Abelian so- 
lutions coincide, one can connect the corresponding pa- 
rameters. The asymptotical value of the potential A$ for 
the Melvin solution is 2/B. In the non- Abelian case we 
have shown that the corresponding asymptotic value is 
—k. This means that 



-Pn 



2/k. 



(165) 



Numerical computations confirm exactly this result, as 
illustrated on Fig. [3J Furthermore, in the Abelian case 
the bifurcation condition QlOSp 



2 = kl 



(166) 



gives again (I165[) . In this limit the following inequality 
holds Eqo > 3 (the strict equality holds for the Melvin 
solution) and thus (3 > — 1. Therefore the parame- 
ter f3 in Eq. (|124H126p takes values within the interval 
[— 1 , 1] , which ensures positivity of the expression under 
the square root. But the signs before the square root are 
different in the Abelian and non-Abelian cases. 

The space generated by the i/3-string can be visualized 
as the surface of a cup. Since < b < 1, the circumfer- 
ence of the circle bounding the disk of the radius r cen- 
tered on the axis and perpendicular to it is growing up 
slower than r. If we embed the surface (r, tp) carrying the 
corresponding metric into the flat space of one dimension 
more, it will look like a cup whose shape looks asymp- 



-X v and the Kasner parameters can be expressed totically as a figure of rotation of the curve Zfi at 



flat 
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Fig. [TDJ Two such surfaces lying at the unit distance 
along the axis will form a thick disk with the thickness 

decreasing asymptotically as z fi at = rj}°\ , see Fig. QTJ 

Solutions corresponding to any set of parameters 
(q, p, k) not belonging to the bifurcation curve are singu- 
lar. Their behavior is similar to that in the flat case, but 
now there is a singularity of the metric, not of the gauge 
field. At the singularity the determinant of the metric 
vanishes, g(r s ) = 0, while the Kretschmann scalar ([70]) 
diverges as (r s — r)~ 4 . The corresponding expansions of 
the metric and the field functions in terms of the devia- 
tion x = r s — r can be found analytically. It turns out 
that the expansion of field functions contains only the 
integer exponents while that of the metric functions - 
only half-integer exponents. In the leading order in x 

P, R -> const, N -> x~ 1/2 , L, K -> x 1/2 . (167) 

An example of such an expansion for k = 1 is given in 
Tab. ED 

At the end of this section we would like to mention 
that the obtained results for the bifurcation conditions 
(|108[ I118j) . as well as the asymptotic expansions (|145I 
I146| I150H154|) does not depend on the assumed discrete 
nature of k. The asymptotic behavior of the system will 
be the same for non- integer values of k. 

VI. CONCLUSION 

In this paper, we have continued investigation of static 
cylindrically symmetric purely magnetic SU(2) EYM 
configurations started in [T^]. We extended analysis of 
bifurcation manifolds, gave a (non-rigorous) proof of ex- 
istence, considered in more detail structure of solution 
space near the high-curvature limit and described in de- 
tail singular solutions. Our set of the integral sum rules 
allows one to relate the asymptotic parameters of the so- 



lution to the boundary data on the polar axis and to sim- 
plify the numerical integration. The values of the Kasner 
exponents can be found analytically in some cases. Oth- 
erwise they are obtained numerically starting with the 
data on the symmetry axis parameterized by the value 
of the magnetic field. 

The geometric structure and the matter distribution 
for soliton solutions obtained suggest an interpretation of 
the solution as describing a pair of straight and circular 
magnetic strings. Asymptotically the metric exhibits the 
boost symmetry in the plane (t, <p) similar to the (t, z) 
boost symmetry of the usual Melvin solution. 

Some generalizations can be suggested, in particular, 
similar structure should exist in the bosonic sector of 
the Freedman-Shwartz SU(2)xSU(2) gauged supergrav- 
ity. An intriguing question is whether the supersym- 
metric solutions may exist. The corresponding Abelian 
U(l)xU(l) configurations were considered in [231 show- 
ing the absence of supersymmetry for the Melvin-type 
configuration. But as was shown in [H, [2§| , spherically 
symmetric static pure magnetic configurations admits su- 
persymmetry of gauged supergravity only in the case of 
truncation to SU(2) xU(l). So supersymmetry for cylin- 
drically symmetric gauge field configurations is still not 
excluded in the non- Abelian sector. 

As other directions we could mention an inclusion of 
the electric component of the YM field, as well as dress- 
ing the solutions with additional structure such as_prop- 
agating gravitational waves [3(| and black holes [3l], [32| , 
similar to that found for the Melvin background. 
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TABLE I: Expansions at the origin in flat space 
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TABLE II: Expansion at the origin in curved space (k — 1) 
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TABLE III: Expansion at the singular point in curved space (k — 1) 
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FIG. 1: Bifurcation of solutions in fiat space. Dashed lines — the two singular solutions corresponding to the parameter 
values p — —0.9076 and p — —0.9078 respectively. Solid lines correspond to the regular solution with the bifurcational value 
p= -0.90777863505.... 
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FIG. 2: Bifurcation curves for k = 1, 2, 3, 4 in the flat space. 
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FIG. 3: Bifurcation curves for k = 1, 2, 3 in the self-gravitating case (solid lines) and in the flat space (dashed lines). 
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FIG. 5: The metric and field functions for the self-gravitating solutions. Dashed lines — the singular solution (q — 0.39, p = 
-0.5), solid lines — the regular solution (q = 0.3862649505 . . . , p = -0.5). 
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FIG. 8: The metric of the regular solution with p = —1.85. Solid lines correspond to the effective exponents of metric functions 
hi — (hxNi)'r where i labels a, b, c and N, L, K respectively. When a% w const, the metric functions Ni ~ r a *. Dashed line 
— the energy density a = T® yj—g. 
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FIG. 9: Continuation of the previous plot Fig. [S] with an enhanced (factor of 100) scale for a. 
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FIG. 10: The geometry of the section z = const, of the solution carrying the metric dl 2 = dr 2 + L 2 (r)dip 2 as embedded into the 
three-dimensional flat space (the desired surface is the figure of rotation of this curve around the vertical axis). The coordinate 
along the horizontal axis is rfi at = yf§ vv , the vertical coordinate is defined by dz 2 lat = dr 2 — dr 2 lat . Three regions correspond 
to the nearly flat metric, and that of the Melvin straight and circular strings respectively. 




FIG. 11: Plot of the metric function K(r) showing an effective longitudinal thickness of the configuration along the z— axis as 
a function of the radial. Three regions are shown which correspond to the nearly flat metric, the straight string core and the 
circular string. 



